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TWO-POTENTIAL APPROACH TO MULTI-DIMENSIONAL 

TUNNELING 



S.A. GURVITZ 

Department of Particle Physics, Weizmann Institute of Science, 
Rehovot 76100, Israel 

We consider tunneling to the continuum in a multi-dimensional potential. It is 
demonstrate that this problem can be treated as two separate problems: a) a 
bound state and b) a non-resonance scattering problem, by a proper splitting 
of the potential into two components. Finally we obtain the resonance energy 
and the partial tunneling widths in terms of the bound and the scattering state 
wave functions. This result can be used in a variety of tunneling problems. As 
an application we consider the ionization of atomic states by an external field. 
We obtain very simple analytical expressions for the tunneling width and for the 
angular distribution of tunneling electrons. It is shown that the angular spread of 
electrons in the final state is determined by the semi-classical traversal tunneling 
time. 
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2 S.A. Gurvitz 



1 Introduction 

This paper is devoted to the memory of Michael Marinov, my friend and 
collaborator. Within the scope of his broad scientific activity Michael paid spe- 
cial attention to different aspects of quantum mechanical tunneling, as for iru 
stance to the coherent and quantumjjiacillation effects in tunneling processespH 
and to the traversal tunneling timeBa In particular Michael was interested in 
a semi-classical description of multi-dimensional tunneling. Our frequent dis- 
cussions of this subject were very illuminating and helpful to me, especially 
in relation with the two-potential approach to multi-dimensional tunneling, 
which I shall present below. 

It is well known that propagation through classically forbidden domains 
has been extensively studied since the early days of quantum mechanics. Yet, 
its exact treatment still remains very complicated and often not practical, in 
particular for the multi-dimensional case. The treatment of the tunneling prob- 
lem can be essentially simplified by reducing it to two separate problems: a 
bound state plus a non-resonance scatterirLg state. This can be done consis- 
tently in the two-potential approach (TPA)EHII This approach provides better 
physical insight than other existing approximation methods, and it is simple 
and very accurate. It was originally derived by us for the one-dimensional 
case. Here we present an extension of TPA to many degrees of freedom and 
its application to ionization of atomic states. 



2 General formalism 

Consider a quantum system in a multi-dimensional space, described by 
the Schrodinger equation (we have adapted units where h = 1) 

i^^\^{t)) = [K+v{xmit)), (1) 

where x — {xi, X2, ■ ■ ■ , Xn}, K = V^/2r7i is the kinetic energy term and V{x) 
is a multi-dimensional potential. The latter contains a barrier that divides the 
entire space into two domains (the "inner" and the "outer" one) , such that the 
classical motion of the system in the inner domain is confined. Yet this does not 
hold for the quantum mechanical motion of this system. It is well known that 
any quantum state with a positive energy, localized initially inside the inner 
domain, tunnels to continuum (the outer domain) through the barrier. If the 
time-dependence of this process is dominated by its exponential component, 
we have a quasi-stationary (resonance) state. Below we describe a simple and 
consistent procedure for treatment of such multi-dimensional quasi-stationary 
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states, which represents a generahzation of the TPA, developed in Rcfs.^-^ 
for one-dimensional case. 

Let us first introduce a hyper-surface S inside the potential barrier, which 
separates the whole space into the inner and the outer regions, as shown 
schematically in Fig. 1. Respectively, the potential V{x) can be represented as 




Figure 1: An example of the separation between the inner and the outer regions for a two- 
dimensional potential. Si and S2 denote the equipotential lines for the energy Eq, Eq. (BI). 
The shaded space correspond to the classically forbidden region. The separation line E is 
taken as a circle. 



a sum of two components, V{x) = U {x) + W{x). The first component U{x) is 
defined in such a way that U{x) — V{x) in the inner region, and U{x) — Uq in 
the outer region, where Uq is the minimal value of V{x) on the hyper-surface, 
Uq =min V{x € E). Respectively, the potential W{x) = in the inner region, 
and W{x) = V{x) — Uq in the outer region. 

We start with the bound state ^o{x) of the "inner" potential, 

(X + C/)|$o> =^o|*o), (2) 

where Eq < Uq. The potential W{x) is switched on a,t t ~ 0. Then the state 
|$o) becomes the wave packet 

/lih 
bi.mUx)e-^'''''j^ (3) 

where bo{t),bk(t) are the probability amplitudes of finding the system in the 
corresponding eigenstates |$o,fc) of the Hamiltonian Hq = K + U. For sim- 
plicity we assume that the potential U{x) contains only one bound state |<i>o) 
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corresponding to the energy Eq, Eq. (|^). The amplitudes bo{t), bk{t) can be 
found from the Schrodinger equation idt\'^) = (iJo + supplemented 
with the initial conditions: 6o(0 = 1j ^fc(O) = 0. 

If the state Y^{t)) is a metastable one, the probability of finding the system 
inside the inner region drops down exponentially: |foo(i)P ~ where T is 

the width of this state. This implies that T is determined by the exponential 
component of the amplitude ho{t)- The latter corresponds to a pole in the 
complex .E-plane of the Laplace transformed amplitude bo{E): 

bo{E)= / e'''Ho{t)dt (4) 
Jo 

The amplitude h^{E) can be found fronuthe Schrodinger equation by using the 
Green's function technique. We obtainJa 

ho{E) = ^ . (5) 

E - E(, - {^o\W\^o) - {^n\WG{E)W\<i>Q) 

where the Green's function G is given by 

G{E)^Go{E){l + WG{E)), (6) 

and 

GniE) = ^—^ , A=|$o)($o|. (7) 

^ ' E + Uo - K -U I u/\ ui w 

Here W{x) = W{x) + Uq is the outer component of the potential V{x). Note 
that Eqs. (|^) represent the standard perturbation theory, except for the 
"renormalized" distortion pijtential, i.e. W{x) is replaced by W{x) in the 
Green's function G, Eq. 

It follows from Eq. (g) that a pole of bo{E) (at E = Er ~ iT/2), that 
determines the energy Er and the width F of a quasi-stationary state, is given 
by the equation 

E = Eo + ($o|W^|$o) + {<^>o\WG{E)W\^o) ■ (8) 

One can easily showi that Eq. (||) determines in fact, the poles in the total 
Green's function G{E) — {E — K — V)^^ at complex E. Since their position 
depends only on the potential V{x)^ the resonance energy Ej. and the width F, 
given by Eq. (H) are independent of the choice of the separation hyper-surface, 
S. 

One can solve Eq. perturbatively by using the standard Born series for 
the Green's function 

G = Go + GoWGo + GoWGoWGo + --- . (9) 
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Substituting Eq. (|) into Eq. (|) we obtain the corresponding series for Er — 
iV /2. Yet, such a series converges very slowly. 

For this reason we proposed in Ref. ^ a different expansion for G, which 
converges much faster than the Born series. This is achieved by expanding G 
in powers of the Greens function , 

G^[E) = i ^, (10) 

^ E-K-W 

corresponding to the outer potential W, instead of Gq, Eqs. (|^), fl), which 
corresponds to the inner potential U . One then obtains from Eqs. ( p^ 

G = G^ + G^[U -Ua)G-G^K{l + WG) (11) 

Iterating Eq. ( pi] ) in powers of G^ and then substituting the result into 
Eq. (|8|) we obtain the desirable perturbative expansion for the resonant energy 
and the width. Using simple physical reasonings and direct evaluation of higher 
order terms for one dimensional case, we demonstrated in Refs.^, ^ that this 
expansion converges very fast, so that G is well approximated by the first term 
Gyf^. In this case Eq. (^) becomes: 

E = Eo + ($o|T4^|$o) + {^o\WG^iE)W\^o), (12) 

Solving Eq. (|lj) for complex E = E,. — iV /2 we find the resonance energy E^ 
and the width F with a high accuracy (of order T /Uq). 

Let us assume that that | A | , F ^ iJg , where ^ = Er — Eq \s the energy 
shift. In this case we can solve Eq. ( |l^ ) iteratively by replacing the argument 
E in the Green's function by Eq: 

E = Eo + (^olT^l^o) + {'i>o\WG^iEo)W\<S?o) ■ (13) 

Equation (|l|) can be significantly simplified. Indeed, the Green's function Gyh 
satisfies the Schrodingcr equation 

Eo-K - W{x) \ G^{E; £c, x') = 5{x ~ x') , (14) 

which can be rewritten as 

W{x)G^{Eq- X, x') = (^Eo -Uo-'^^ G^{Eo; x, x') - S{x - x') . (15) 

Substituting Eq. (^5|) into Eq. (O) we find that the corresponding contribution 
from the kinetic energy term (V /2to) can be transformed as 

Mx)V^G^{Eo; X, x') = V{<^oVG^) - V{G^V'i>o) + G^^ V^<i>o (16) 
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In addition, the potential W{x) is zero in the inner region, so that the integra- 
tion over X, x' in the matrix elements (|l3|) takes place only in the outer region, 
TZout- Then using Eq. (|l6|) and the Gauss theorem we find 



+ / G^V'^^odx (17) 

Since the wave function of the (inner) bound state in the outer region satisfies 
the Schrodinger equation with the potential Uq, the last term in Eq. ( p7| ) can 
be rewritten as 

/ G^V^^odx = 2m([/o - Eo) [ G^^^odx , (18) 
Using Eqs. (|l5|), (|l^, (^ we then obtain 

j M=^)W{x)G^{Eo;x,x')dx = ^ J [Gi^,V„$o-$oV„G^]da-$o(a;') 

(19) 

Substituting Eq. (p3) into Eq. (p^ we find that the latter can be rewritten as 



E^Eo-:^ [ Mx) V„ Gyy{Eo;x, x')W{x')<Po{x')dada' , (20) 
■^''^ JxeT. 

where V means the gradient on the right minus the gradient on the left. 

The x' integration can be carried out in the same way by using Eqs. (|l5|)- 
(|l8|). Again, the V^$o(3;') term cancels as well as the (5-function contribu- 
tion. Finally we obtain that the resonance energy and the width are expressed 
through a double integral over the hyper-surface E 

^^^^-TT^I Mx)Vr,G^{EQ-x,x')Vr,' '^o{x')dada' . (21) 

The resonance width T can be easily obtained from Eq. ( ^l| ) by taking the 
spectral representation for the Green's function, 

„ _ /■ |yfc)(yfc| dk 
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where Ek = k /2m and (pk are the eigenstates of the "outer" Hamihonian 
correspondmg to non- resonance scattermg states. 



[K + W{x)]ipkix) = Ekipk{x) 



(23) 



Using Eqs. (^l|), ( p^ ) one obtains the total width as an integral over the partial 
width Ffc 

dki ■ ■ ■ dkn-i 



r= r 



where ko = \/2mEo and 



(27r)"-i 



(24) 



|fc|=fco 



Tk = 



Amk„ 



(25) 



\k\=ko 



Equations. (|24| ) and (g5|) represent our final result for the total and partial 
tunneling widths of the resonance state with energy Er — Eq, Eq. (||) The 
latter is the energy of a pure bound state 4>o in the inner potential U{x). Both 
Tk and Eq are independent of the hyper-surface S (up to the terms of the order 
of T/Uq), providing that E is taken far away from the equipotential surfaces 
V{x) ^ Eo B 

Equations. (||), (|2^), ( p5| ) considerably simplify the multi-dimensional tun- 
neling problem, by reducing it to the bound state problem. The latter can 
be treated with standard techniques. Some examples of applications of our 
method are given below. 



3 Ionization of bound state in semi-classical approximation 

Consider a three-dimensional potential u(r), which contains a bound state 
(/)o(r) with the energy eg < 0. Let us superimpose an additional potential v{x), 
depending on the coordinate x. We assume that v{x) ionizes the bound state 
(/)o, i.e. u{r) + v{x) < eq for x ^ oo. We also assume that u{r) and v{x) do 
not overlap, such that u{r) « for x > xi and v{x) ~ for x < X2, where 
xi < X2- Now taking the separation hyper-surface S as a plane, x = xq where 
xi < xq < X2, we can identify u{r) and v{x) with the inner and the outer 
potentials, U and W, Eqs. (^), (|2^). Then the bound state wave function 
(j)Q{r) coincides with the inner wave function $o; Eq. (||). 

Providing that v(x) is a function of x only, the outer wave function, 
Eq. (^) can be written as a product fkir) = eyi-p{iqb)xk{x), with r = {x,b), 



° Equation (|25) resembles the Bardeen formulaO for energy splitting of two degenerate levels 



in a double-we 



1 potential. 



8 S.A. Gurvitz 



k = {k,q) and = 2m(uo + eo)j where uq = —v{x oo). As a result, 
Eq. ( p5| ) for the partial width Tq reads 

1 



Arak 



d d 
Mx:b)—x-k{x) - x-k{x)-Q^(t>o{x,b) 



(26) 
I, is 

-(2^ (2^) 

Since Xhi^) ^ non-resonant wave function, it falls down exponentially 
inside the barrier and can be evaluated by using the standard semi-classical 
approximationia One can write 



where k — ^j2m(uQ + eq) ^ 1^- Respectively, the total width, Eq. 
given by 



\Pqix)\ 



ff \pq{x')\dx' 



(28) 



Here Pq{x) = •\/2m[eo — v{x)\ — q^^ and Xf\s the outer classical turning point, 
Pq{xf) = 0. 

The inner wave-function describes the bound state and therefore it de- 
creases exponentially under the barrier, as well 



00 (a;, 6) 



* ' \p(x',b)\dx' 



(29) 



Here \p{x, b)\ = ^/2m\eQ — u{x, b)], and Xi is the inner classical turning point. 
It follows from Eq. ( P9|) that (j)o{x, b) decreases very quickly with b so that only 
small b contribute in the integral (|2g). Thus, we can approximate the inner 
wave function and its derivative as 



Xi - \p{x',0)\dx'-\pix.,0)\b^/2x 

(j)o{x, b) ~ — (po(Xi,0) e •'^i 

X 



(30) 



d(t>Q{x,b) 
dx 



-\p{x,Q)\ Mx,b) 



(31) 



Substituting Eqs. (p8|), (^0|), ( ^ into Eq. ( p6| ) we can perform the b- integration 
by using 



' (l)o{xo,b)d^b 



27^a;i , , ^, - f"" |p(a;',0)|da:'-a:oqV2|25(a;o,0)| 



|p(a;o,0)| 



(32) 
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Since xq 3> Xi the exponent in Eq. ( p2[ ) can be approximated as 

Now using Eqs. (|3^), (|33|) , and dxkl^^^ — \Pq{'^)\Xk^ obtain the foUow- 
ing formula for the partial width, 

m(2m\eo\) ' 

where 

Pg(a;) = v'2TO[u(a:, 0) + v{x) - eo] + , (35) 

is the absolute value of the imaginary momentum in the direction of the ex- 
ternal field. 



4 Partial width and the tunneling time 

In contrast with the one-dimensional tunneling, the final state in a multi- 
dimensional tunneling is highly degenerate. Hence, a tunneling particle can 
be found in different final states with the same total energy. In our case these 
are all the states with different lateral momenta q, leading to the angular 
distribution of a tunneling particle. Using our result for the partial width, 
Eq. (p4[), one easily finds that this distribution is of a Gaussian shape. Indeed, 
by expanding the exponent of Eq. (BJ) in powers of we obtain 



Pq{x)dx ~ S{EQ,x^,Xf) + — I -r^, (36) 

2"! J,^. p{x) 

where p{x) — Pq=o{x) and 

S{eo,x^,Xf) = / p{x)dx (37) 



is the action in the classically forbidden region. Then Eq. (34) reads 
fn{2m\eo\ )'-''' 

where 

\p{x)\ 
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is the Biittiker-Landauer tunneling tinic0 The latter is assumed to be an 
amount of time, which a particle spends traversing a tunnel barrier. It ap- 
pears, however, that such a quantity is not uniquely defined in one-dimensional 
tunnelingEila Yet, the situation is different for the multi-dimensional case. Our 
result, Eq. (^8|), shows that the Biittiker-Landauer traversal tunneling time r 
is related to a well defined quantity, namely to the lateral spread of a tunnel- 
ing particle (ro).u Indeed, the latter corresponds to the inverse average lateral 
momentum, ro = 1/ < (? >, where < q >= ^fmjr as obtained from Eq. (|3^). 



The Gaussian shape of angular distribution of tunneling particles is a gen- 
eral phenomenon, which can be understood in the following waj'u It is well 
known that the most probable tunneling path corresponds to the minimal ac- 
tion. In our case this path is along the direction of the external field. Therefore 
any momentum q, normal to this direction raises the action by a term propor- 
tional to the corresponding energy /2m. This leads to the Gaussian fall-off 
of the related probability distribution. 



5 Total width and angular distribution of ionized particles 



Substituting Eq. (38) into Eq. (27) and performing the g-integration, we 
obtain the following expression for the total width 



2S{Eo,Xi,x f) 



(2m|eo|)i/2 

can be rewritten in a more compact form as 
47r 



Using this result, Eq. ( 

m 

Let us consider two exemplifying applications of Eqs. ( ]4C| ) 
(a) Hydrogen atom 



(40) 



(41) 



Consider the Hydrogen atom in the ground state under the electric field £. 
We use the atomic units in which m — 1, the ground state energy eo = —1/2 
and the ground state wave function 4>o{r) = exp(— r)/Y^. In these units the 
inner and the outer potentials are: u{r) — ~l/r and v{x) — —Ex. Then the 
action S, Eq. ( ^ reads 

/■^/ j 2 f'^f . dr 

/ \ l 2£xdx^ / VI - 2£xdx - / — , (42) 

Jx, V X J^. J^. xVl - 2Ex 
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Although the main contribution to the action is coming from large x, we have 
keep in (42) also the next term of expansion (see Ref . p^) . Performing the 
integration in Eq. (^) we obtain for the action 

where we assumed that Exi <C 1. Respectively for the "tunneling" one obtains 
T = \/£. Using these results we easily find froroJilq. ( ^o|) the well-known result 
for the ionization width of the Hydrogen atomO 

r = ^ e-2/3f . (44) 

Respectively, the angular distribution of ionized electrons, Eq. (^l]), is 

r, = fre-^/^. (45) 

(h) Potential well with short range forces 

The bound s-state wave function of a short range potential at distances 
larger then the radius of the well a has the following form 

^o(r) = ^e— , (46) 
r 

where k = (2|eo|)"'^/^ ajid ^ is a dimensionless constant depending on the 
specific form of the wellt2l. For instance, if a ^ 1/k, then A — (27r)^^/^. 

Similar to the previous case we use the same units but now with k — 1 . The 
action and the tunneling time are evaluated straightforwardly from Eqs. (^7|), 
( ^9| ) taking into account that the initial turning points are Xi = a — *■ 0. One 
finds 

l/2£ 1 

Vl - 2£xdx = — (47) 



3f 

1/2^ dx 1 



(48) 



y/l - 2£x E 
Using xl4>Q{xi, 0) — one obtains from Eq. ( pO| ) 

r = tttI^ f e-2/3£ _ (49) 
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This coincides with the resuh obtained by Demkov and DrukarevPlEj using a 
more comphcated technique. 

The angular distribution ol ionized particles, Eq. (^l|), is 

= 4^2^2g-2/3£e-'j'/^ . (50) 

The same result was obtained earlier in Ref. IJ. 

6 Summary 

In this paper we extended the two-potential approach to quantum tunnel- 
ing to a multi-dimensional case. We found simple expressions for the energy 
and the partial widths of a quasi-stationary state, initially localized inside a 
multi-dimensional potential. The partial width, which determines the distribu- 
tion of tunneling particles in the final state, is given by an integral of the bound 
and non-resonant scattering wave functions over the hyper-surface taken inside 
the potential barrier. 

We applied our results for the ionization of atomic states by an external 
field depending only on one of the coordinates. In this case we obtained very 
simple semi-classical expressions for the total ionization width and the angular 
distribution of tunneling electrons, valid for any initial state. We demonstrated 
that in the case of the Hydrogen atom and also for a short range potential our 
general expression for the total width produces the results, existing in the 
literature, but obtained with a more complicated calculations. 

We also demonstrated that angular distribution of tunneling electrons in 
the final state falls off as a Gaussian function of the lateral momentum. The 
exponential is determined by the traversal semi-classical tunneling time. 
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